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Abstract. Stochastic resonance with white noise has been well established as a potential signal amplifica- 
tion mechanism in nanomechanical two-state systems. While white noise represents the archetypal stimulus 
for stochastic resonance, typical operating environments for nanomechanical devices often contain different 
classes of noise, particularly colored noise with a 1// spectrum. As a result, improved understanding of 
the effects of noise color will be helpful in maximizing device performance. Here we report measurements 
of stochastic resonance in a silicon nanomechanical resonator using 1// noise and Ornstein-Uhlenbeck 
noise types. Power spectral densities and residence time distributions provide insight into asymmetry of 
the bistable amplitude states, and the data sets suggest that 1//" noise spectra with increasing noise 
color (i.e. a) may lead to increasing asymmetry in the system, reducing the achievable amplification. 
Furthermore, we explore the effects of correlation time r on stochastic resonance with the use of exponen- 
tially correlated noise. We find monotonic suppression of the spectral amplification as the correlation time 
increases. 

PACS. 85.85.+j MEMS/NEMS - 05.40.-a Fluctuation phenomena, noise, and Brownian motion - 05.45.-a 
Nonlinear dynamics and chaos 



1 Introduction 

Although it is not surprising given the wealth of mechani- 
cal, electronic, biological and other systems which display 
stochastic resonance (SR) yij, the observation of noise- 
enhanced switching in nanoelectromechanical systems 
(NEMS) holds promise for the possible exploitation of 
noise to improve device performance [2 . In addition, the 
discovery of SR in NEMS has made available yet another 
experimental platform for the study of basic questions 
in the field, such as the effect of noise color. Combining 
these motivations, we discuss measurements of SR in a sil- 
icon nanomechanical resonator using various colored noise 
spectra. 

Improving the model for SR in this complex applica- 
tion warrants additional attention as NEMS continue to 
emerge as a viable technology. In particular, the bistable 
amplitude states of a nonlinear nanomechanical resonator 
may see use as two-state devices such as memory elements 
and switches f3','T. A particle subject to white noise in 
a periodically-modulated double-well potential remains a 
powerfully intuitive model for stochastic resonance, cap- 
turing the basic dynamics in a diverse array of systems. 
Still, much work has been devoted to generalizing this pic- 
ture. Asymmetric potentials [5/6,T,'5], spatially extended 
systems [9], and "non-conventional" stochastic resonance 
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between coexisting periodic attractors [TOl have been among 
the many extensions to the original idea, with the lat- 
ter two offshoots both applying to nanomechanical SR. 
Such derivations continue to be important in understand- 
ing stochastic resonance in actual practice, where the ex- 
perimental conditions are rarely straightforward. 

Yet another important and ongoing consideration is in 
understanding the effect of noise color on stochastic reso- 
nance [TT | ll2 p i3 } [T4 t [T5 l ll6 p i7; . In particular, the ubiquity 
of 1// noise makes it a prominent factor in a wide range 
of settings [TF , including in the active electronic elements 
used to drive nanomechanical resonators |19j . This noise 
class has been used to induce SR in electronic [13], neu- 
rological |15j . and most recently, nanoelectromechanical 
systems [T7]. Although the suppression of SR persists as 
a common theme in these studies, to the authors knowl- 
edge, an accepted physical picture for the effect of 1//" 
noise exponent is lacking. Asymmetry has been postulated 
as a consequence of 1// noise, and evidence in this study 
further reinforces that idea. 

Extending the study to another class of colored noise, 
we also observe SR with the use of exponentially corre- 
lated (Ornstein-Uhlenbeck) noise. We find that the am- 
plification of the signal declines with increasing correla- 
tion strength of the noise, in agreement with theoretical 
predictions and previous experimental results jll i,.12, 14 . 

Previous studies of SR in nanomechanical resonators 
have focused mainly on traditional measures, such as signal- 
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to-noise ratio (SNR) and spectral amplification [2l[T7ll20|. 
Given the technological potential of these systems, en- 
hancements seen in SNR and signal amplification provide 
obvious implications for the utility of SR in an applied 
setting. However, the wealth of analyses that the field has 
grown to embrace, including of power spectral densities 
and residence time distributions, afford further insights. 
We use these approaches here to discern the role of asym- 
metry on system dynamics. 



2 Experimental Methods 

2.1 Device Fabrication and Characterization 

Using e-beam lithography, resonators are patterned on the 
single crystal silicon layer of a silicon-on-insulator (SOI) 
wafer. The length L and width w, 20 /im and 300 nm, 
respectively, place the in-plane resonance frequency of the 
doubly-clamped beam in the MHz range. A gap d of 250 
nm at equilibrium separates the resonator from adjacent 
parallel electrodes used for actuation and detection, and 
the device layer thickness is 500 nm. Metal contacts are 
deposited via thermal evaporation, giving the resonator a 
thickness t of 550 nm. Finally, reactive ion etching (RIE) 
and hydrofluoric acid etching are used to define and sus- 
pend the silicon structure. 

To actuate the beam, we utilize the standard electro- 
static technique at room temperature ([3] and references 
therein). With a DC bias Vb charging the resonator, a 
high frequency excitation signal vu(t) applied to the drive 
electrode creates a capacitive force. The system can be 
modeled as a Duffing oscillator [21] : 



X + '-fx + uJqX + k^x^ — 



1 dC 
2 



VBVD{t) 



(1) 



where x is the effective resonator displacement, 7 is a dis- 
sipation coefficient, ujq is the resonant angular frequency 
of oscillation, and is the non- linear coefficient. Here, 
the resonator and drive electrode may be accurately ap- 
proximated as parallel plates, making the capacitance be- 
tween them C = eoLt/ (d-x). In-plane oscillation of the bi- 
ased resonator induces a current at the opposite electrode, 
which we amplify and measure with a network analyzer. 
The frequency response exhibits the expected Lorentzian 
line shape for small displacements (Fig. [l|a)). As drive 
power is increased, nonlinearity leads to spring harden- 
ing, and two stable amplitude states emerge with the on- 
set of hysteresis. These bistable amplitude states provide 
a two-state system for the study of SR. 



2.2 Modulation and Noise 

Whereas previous studies of SR in nanomechanical sys- 
tems have utilized an additional periodic signal for switch- 
ing "2 ,20], we modulate the potential by means of a phase 
deviation (po in the drive signal itself. Recent experiments 
have shown that an abrupt phase shift can be used to in- 
duce a switch between the bistable amplitude states of a 
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Fig. 1. a) Resonator response vs. frequency at various drive 
levels. Beginning as a Lorentzian with resonance frequency /o 
~ 3 MHz and quality factor Q « 50, the resonance bends to 
higher frequency with increasing drive voltage. Eventually, up- 
ward (empty squares) and downward (filled circles) frequency 
sweeps at a single drive voltage exhibit hysteresis, revealing two 
stable amplitude states, b) Micrograph of the resonator with 
SR measurement setup. Separate sources provide the modu- 
lated drive signal f£)(t) and the noise spectrum ^{t), while a 
network analyzer operating in continuous wave mode monitors 
the response amplitude at the drive frequency lo. 



nanomechanical resonator [4]. Here, we drive the system 
in the bistable region and modulate the phase of the drive 
with a square wave: 



U£)(i) = VBCOs{ujt 



(2) 



where 0{f2t) represents a square wave of period ^, alter- 
nating between -1-1 and -1 each successive half period. The 
drive frequency uj /2t: is a frequency in the bistable region, 
and J7 << aj/27r (J7 = 50 Hz for all results reported here). 
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Expanding this expression, we see that the drive consists 
of two parts: 

voit) — VD[cos{'^)cos{ujt)+0{f2t)sin{^)sin{u}t)]. (3) 

The first term in the brackets provides simple periodic 
forcing. The second term, which changes phase by 180 de- 
grees each half period of the square wave, induces switch- 
ing of the resonator amplitude in synchronization with 
0{f2t). The magnitude of ipo dictates the amplitudes of 
these two terms, and for small angles, the amplitude of 
the modulation signal is proportional to ipo. The exact 
effect of this modulation depends upon the specific drive 
conditions of the resonator, but in general, phase devia- 
tions close to TT radians induce a switch. As the magnitude 
of the phase deviation is dropped, the switching fidelity 
drops as well, and below a certain value of ipo, the mod- 
ulation becomes sub-threshold and the resonator remains 
in the given amplitude state (see, for example. Figure 3 in 
In these experiments, we set the periodic modulation 
just below this threshold and recover coherent switching 
via the addition of relatively small noise intensities to the 
system. This places the measurements in the weak noise 
limit, outside the applicability of linear response approxi- 
mations 

Colored noise spectra are digitally generated in Mat- 
Lab and reproduced by a second signal generator. The sig- 
nal generator implements I/Q modulation of a carrier sig- 
nal, producing a voltage ^{t) — I{t)cos{uJct)+Q{t)sin{ujct) , 
where I(t) and Q{t) are the digitally generated noise se- 
quences, and ujc is the angular frequency of the carrier 
signal. We generate 1//" noise by creating a sequence of 
63,000 points in frequency space with an amplitude of the 
form 1// 2" and random phase (taken from a uniform dis- 
tribution [— TT, tt]). Taking the inverse Fourier transform of 
this set, we obtain the desired noise sequence in the time 
domain. In the case of Ornstein-Uhlenbeck noise, the noise 
sequences are generated by solving the equation: 

iou{t) = —iou{t) + —^n,h{t) (4) 
T T 

where r is the correlation time, D is the noise intensity 
and £,wh{t) is a sequence of 63,000 points of Gaussian white 
noise generated in MatLab using the randn{) function. 

2.3 Measurements 

Bistability of the resonator extends over a relatively wide 
range of operable frequencies and drive voltages. Further- 
more, the drive parameters influence the relative stability 
of the two amplitude states [5T], affecting the symmetry 
of the potential describing the system. As such, care is 
taken prior to measurement to symmetrize the states as 
closely as possible. In the absence of noise and with the 
phase modulation inducing infrequent switching, we ad- 
just the high frequency drive voltage until neither ampli- 
tude state is strongly preferred. The phase modulation is 



then dropped to a sub-threshold level, and noise is added 
within a 50 kHz band encompassing the full width at half 
maximum (FWHM) of the resonance. In the case of 1// 
noise, we choose an I/Q carrier frequency at the low fre- 
quency side of the beam resonance, so that noise power 
spectral density (PSD) decays across the bistability region 
in frequency space. By contrast, we center the Lorentzian 
line shapes of exponentially correlated noise within the 
bistability region itself. 

Figure [ijb) displays the measurement setup. For each 
value of noise intensity, a network analyzer monitors the 
response at the drive frequency uj for 500 periods of the 
modulation. We measure SR using 1//" noise with a 
ranging from to 2 in steps of 0.5. Using exponentially 
correlated noise, we measure SR with noise correlation 
times spanning nearly three orders of magnitude, from r 
w 0.01i7-i to T w 10/2-1. 



3 Results 

3.1 1/f Noise 

Trends in the SNR using 1/f noise are reported elsewhere 
[T7]. Specifically, results show that increasing the noise 
exponent suppresses the resonance in SNR, shifting it to 
higher noise powers and reducing the peak values reached. 
We note that interpretations of the underlying cause of 
these trends may vary, with an increasing effective po- 
tential |16l and an increasing system asymmetry [5] both 
serving as potential explanations. Here we focus on al- 
ternative analyses that may help to illuminate the effects 
governing SR with 1/f noise. In these analyses, we draw 
heavily upon theoretical works which treat the canonical 
Duffing potential. Although the nonlinear behavior of the 
resonator is properly described by Eq. [T] we note that 
application of the double-well model here is not strictly 
correct. The system under study is dynamic, with ampli- 
tude as the proper state variable (not position). Although 
the most accurate treatment would involve coexisting pe- 
riodic attractors |10| . we focus on the canonical model in 
the interest of gaining some physical intuition for the effect 
of noise color on the system. 

To begin, we apply a two-state filter to the recorded 
data and take the Fourier transform to examine the spec- 
tral response. Figure 2 displays the result of a typical mea- 
surement at low noise intensity. Several features present 
in this example arise consistently throughout the mea- 
surements taken with 1/f noise. Most prominently, the 
first few odd harmonics of Q are present, decaying ap- 
proximately as l/'n?, as expected. Also evident are even 
harmonics of the modulation frequency, which show little 
to no dependency on n. Finally, we note the presence of 
broad holes in the vicinity of even harmonics. The latter 
two observations are consistent with expectations for the 
weak noise limit [22]. We note also that the presence of 
holes and/or peaks at even multiples of have been a 
prominent feature of SR studies featuring asymmetric po- 
tentials [HllllT]. At higher noise powers, the holes "wash 
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Fig. 2. An example of power spectral density for a single mea- 
surement. Here, 1//" noise is used with a = 1 and Pn / Pd ~ 
0.01. Squares denote the odd harmonics and circles the even 
harmonics of the modulation frequency. Provided linear fits 
show that odd harmonics decay approximately as while 
even harmonics show no notable dependency with n. 




Fig. 3. Dependency of the harmonics Q„ of modulation fre- 
quency on a for a single noise power, Pn /Pd ~ 0.0065. Here, 
we omit q = 2, as this noise spectrum proves insufficient to 
recover switching at this noise power. 
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Fig. 4. Example residence time distributions for the low and 
high amplitude states, denoted 1 and 2, respectively, at noise 
power Pn / Pd ~ 0.005. Plots display data for: (a) a = 0, (b) 
a = 0.5 and (c) a = 1. 



out" , disappearing as the overall noise level in the spectral 
response rises. 

In general, limitations on sampling rate make quanti- 
tative analysis of the harmonics Qn and the holes difficult 
and error prone; however, we note that at low noise pow- 
ers, extrapolation of the linear fits shown in Fig.[2]provides 



some evidence that the magnitudes of odd and even har- 
monics will cross sooner with higher a. Theoretical treat- 
ments predict oscillation in the relative magnitudes of odd 
and even Qn [7 . As asymmetry increases, the magnitudes 
of the even harmonics surpass those of the odd harmonics 
at lower frequencies, the effect we see here with increasing 
a. In addition, at low noise powers, odd harmonics decay 
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monotonically with the noise exponent (Fig. |3]). By con- 
trast, even harmonics show a very shght maximum as a 
function of a, an effect seen with increasing system asym- 
metry in both numerical simulations [6 and theoretical 
calculations [J. 

Residence time distributions help to provide a more 
complete picture of the system dynamics (Fig. [4]). At low 
noise powers, distributions consistently exhibit preferen- 
tial occupation of one of the states, specifically the low 
amplitude one. Upon reaching the high amplitude state, 
the resonator frequently makes a correlated switch, re- 
turning to the low amplitude state after half a period of 
the modulation. By contrast, occupation of the low am- 
plitude state may persist for multiple periods. In the plot 
shown for a — I, for example, the resonator often remains 
in the low state for durations longer than those pictured, 
with the longest being for more than thirty periods of the 
modulation. Some asymmetry is also evident in measure- 
ments made with white noise, despite previously described 



efforts to symmetrize the potential (Sec. 2.3) 



To quantify any evolving asymmetry with noise color, 
we plot the difference AT between the average residence 
time of the two states vs. input noise power (Fig. [s]). The 
effects of the asymmetry are most pronounced at low noise 
intensities, resulting in significant deviations between the 
observed average occupation times. In all cases, as noise 
power is increased and uncorrelated switching begins to 
dominate, the effects of asymmetry are mitigated, and the 
difference AT decays toward zero. Of note, however, is the 
shift in the AT vs. noise power data as we increase the 
noise exponent a. Theoretical studies of switching time 
distributions predict the following form for AT with an 
asymmetry added to the Duffing potential [7]: 



AT (X 



ce ^« 



3/2 



(5) 



N 



where Pjv is the noise power, AU is the potential bar- 
rier separating the states, and c is the linear coefficient in 
the system potential, essentially determining the degree 
of asymmetry present. Here we express the noise power 
in units of the phase- modulated drive power Pq. Since 
increasing either or both of c and AU could create the 
trend, discerning the cause of the observed shift is diffi- 
cult. In two-parameter least squares fits (not pictured), 
the value of c increases steadily with a while AU exhibits 
no clear trend; however, in some cases, initial parameter 
choices affect the final values. 

For comparison, we also attempt two separate sets of 
one-parameter fits to Equation [5] With AU fixed at its 
value from the white noise (a = 0) two-parameter fit. Fig- 
ure p[a) displays regressions with c as the fitting param- 
eterTlikewise, Figure |5jb) shows fits with AU allowed to 
vary and c fixed at its value from the white noise data. 
Clearly the former set better captures the trend seen with 
increasing a. In essence, a steady increase in asymmetry 
is necessary to account for the severe shift observed. We 
conclude that increasing noise color intensifies the system 
asymmetry, though we cannot rule out the possibility that 
the effective potential barrier deepens as well. 




Fig. 5. Difference AT between the average residence time of 
the two amplitude states vs. input noise power. The difference 
is expressed in units of the modulation period, 1/ Q. The shapes 
denote the type of 1//" noise used: a = (empty squares), a 
— 0.5 (empty circles), a — 1 (empty triangles), q = 1.5 (filled 
squares), and a — 2 (filled circles). In both (a) and (b), the 
line for white noise represents the same two parameter fit to 
Eq. [5] Otherwise, solid lines in (a) and (b) depict fits with c 
and AU allowed to vary, respectively, with the other variable 
held constant. 



3.2 Ornstein-Uhlenbeck Noise 



Extending our measurements to another class of colored 
noise, we also examine the effect of correlation time on 
stochastic resonance with the use of exponentially corre- 
lated noise. Noise power is defined as the power con- 
tained within the FWHM of the resonance, and it is again 
expressed in units of the phase-modulated drive power Pd ■ 
Spectral amplification, defined as 5'(J7)|p„/S'(/2)|p„=o, 
where S{f2) is the PSD at the modulation frequency, is 
monotonically suppressed as the correlation time of the 
noise spectrum becomes longer (Fig |6| . Peak amplifica- 
tion decreases and shifts to higher noise powers as the 
noise color increases. 
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Since it is more easily approached mathematically than 
1// noise, Ornstein-Uhlenbeck noise is more prevalent in 
theoretical treatments of SR. In particular, it has been 
explored via the Unified Colored Noise Approximation 
(UCNA) |T2]. Our measurements in the nonlinear regime 
are beyond the applicability of such treatments, but as a 
means of emphasizing the trends, we fit with two parame- 
ters to the basic form predicted for spectral amplification 
rj: 



a) 



v = i 



A 



Pn/Pl 



1 



(6) 



Solid lines in Figure [6j^ a) depict these fits, from which the 
peak value of rj and the corresponding noise power at the 
peak are extracted and plotted (Fig. [6]jb)). Over almost 
three orders of magnitude in r we observe a very slight 
shift in peak location, a notable contrast with the nearly 
linear shift predicted by linear response theory |12j and 
with the much more significant shift seen when increasing 
1// noise color [17j. 



4 Discussion 

Stochastic resonance using 1// noise has been studied ex- 
perimentally in electronic |13| and biological systems |15j, 
and it has been the focus of some theoretical work [16J. Of 
particular note are parallels to measurements made with 
a Schmitt trigger, in which even harmonics and holes were 
evident in the PSD, and the 1// noise was interpreted as 
creating a slowly- varying asymmetry in the system poten- 
tial |13j . Here we extend the approach by varying the value 
of the noise exponent a and find evidence that increasing 
the noise color exacerbates the asymmetry, an intuitive 
result. As a consequence of the increasing asymmetry, the 
resonance in signal-to-noise ratio should shrink and shift 
to higher noise powers [5], as observed elsewhere [Tf] . 

Using exponentially correlated noise, the suppression 
of stochastic resonance we observe with increasing correla- 
tion time agrees well with theoretical predictions |12j , ana- 
log simulations |rij and experimental results obtained in 
a tunnel diode [Hj . Measurement in the nonlinear regime 
separates these results slightly from those of other inves- 
tigations, however, the conclusion that the use of colored 
noise reduces the effectiveness of SR is a common one. 
To optimize the performance of the nanomechanical res- 
onator as a two-state system, white noise proves superior 
to both 1// and Ornstein-Uhlenbeck noise types, inducing 
greater peak amplification at smaller noise powers. 

Analyzing the results with Ornstein-Uhlenbeck noise 
for trends in asymmetry, we find no evidence of a connec- 
tion similar to that seen with 1// noise, a difference that 
may stem from our placement of the noise spectra in fre- 
quency space (Sec. 2.3 1. Even harmonics and holes in the 
PSD again persist for small noise intensities and disappear 
as the weak noise limit is eclipsed; however, plots of AT 
vs. noise power exhibit no measurable shift with changing 
correlation time. This result, combined with the lack of a 
significant shift of the spectral amplification resonance to 
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Fig. 6. a) Spectral amplification vs. noise power using 
Ornstein-Uhlenbeck noise with three exemplary values of cor- 
relation time T. Additional values of r are omitted for clarity. 
Solid lines depict fits to Eq. |6] Correlation times are measured 
using a spectrum analyzer and are presented here as a unit- 
less number scaled by the modulation frequency (i.e. I/O — 
1). b) Maximum spectral amplification and corresponding noise 
power necessary to reach that maximum, respectively, vs. noise 
correlation time. Values are taken from the fits described in 
(a). Solid lines are included as a guide to the eye and are not 
intended as fits to any theoretical form. 



higher noise powers, suggests no strong effect of the noise 
correlation time on potential depth or asymmetry for this 
system. 



5 Conclusion 

Stochastic resonance is measured in a bistable silicon nano- 
mechanical resonator using a wide range of colored noise 
spectra. Evidence suggests that escalation of 1//" noise 
color increases an inherent potential asymmetry in the 
nanomechanical two-state system, explaining the suppres- 
sion of signal-to-noise ratio seen with this noise class. In- 
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creasing correlation time with Ornstein-Uhlenbcc;k noise 
monotonically suppresses the stochastic resonance, though 
the shift observed in resonance over three orders of mag- 
nitude in T is substantially less than that observed for 
1// noise when changing a from to 2. In either case, 
increasing the color of the noise spectrum weakens the 
effectiveness of stochastic resonance, making white noise 
the best option for optimizing device performance. 

This work is supported by NSF (DMR-0449670). 
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